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Abstract. This overview article concerns the notion of fractional smoothness of 
random variables of the form giXr), where X — (-Xt)ts[o,T] is a certain diffusion 
process. We review the connection to the real interpolation theory, give examples and 
applications of this concept. The applications in stochastic finance mainly concern 
the analysis of discrete time hedging errors. We close the review by indicating some 
further developments. 



1.1 Introduction 

From the practitioners one learns that hedging an option which payoff is 
discontinuous is more difficult compared to the case of smooth payoffs: this 
feature appears for instance for digital options or barrier options (we refer 
the reader to [Tal97] among others). On the one hand, for such options the 
number of assets (i.e. the delta) to incorporate in the hedging portfolio is 
unbounded, and it may become larger and larger when one gets close to the 
singularity (i.e. the maturity and the strike for digital options, or the trigger 
level for barrier options). On the other hand, the numerical estimation of this 
delta becomes less and less accurate, leading to global stability issues. These 
heuristic observations are the starting point for deeper mathematical investi- 
gations about the concept of irregular payoffs, in order to formalize it and to 
quantify the payoff irregularity (with the notion of fractional smoothness). In 
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the current contribution, we aim to give an overview of this concept and some 
applications in stochastic finance. Actually, the applications go beyond the 
financial framework and more generally, they concern the theory of stochastic 
differential equations and their approximations. 

The discrete time hedging error as an important application. 

Since most of the results presented here are applied to the aforementioned 
example of hedging possibly irregular options, we start with a brief presenta- 
tion of this problem, in order to emphasize the issues to handle and to raise 
some natural questions. Take for instance an European-style option exercised 
at maturity T > 0, with a payoff of the form H(St) where S t ■— [Sj, ■ ■ ■ , Sf] 
denotes the price of a d-dimensional underlying asset at time < t < T. 
Sometimes, we will use the notation X\ = log(S'j) for the log-asset, and 
g(x) = h(e Xl , • • • , e Xd ) for the payoff in the logarithmic variables. In what 
follows, we assume a Markovian dynamics without jumps for the asset (solu- 
tion to a SDE defined below), we suppose that the interest rate is equal to 
(to simplify the presentation) and that the market is complete (for details 
about this standard framework, see [KS98]). Thus, under some regularity as- 
sumptions, the payoff h{Sx) can be replicated perfectly by a continuous time 
strategy, where Sf = V x H(t, S t ) defines the vector of number of assets to hold 
at time t. Here, H is the fair price of the option, that is 



where Q is the (unique) risk-neutral measure. In practice, only discrete-time 
hedging is possible at some times t = (U)f =0 with = to < ti < • • • < t n -i < 
t n = T. Thus, at time t € [0, T] the option seller is left with the tracking error 



C t (h(S T ),r) = H(t, S t ) - H(0, So) - £ V X H(U, S u ) ■ (S ti+lA t - S ttM ) 



with <j>(s) = ti when ti < s < U + \. We expect the tracking error (1.1) to 
converge to as the number n of re-balancing dates goes to infinity. With 
the above formulation (1.1), the tracking error is naturally associated to the 
problem of approximation of a stochastic integral using piece-wise constant 
integrated processes. But the delta process {V x H{s 1 £s))se[o,T) may exhibit 
very different behaviors from payoff to payoff: if the payoff is smooth enough, 
then the delta might be bounded as time goes to maturity, while an irregular 
payoff usually yields an exploding delta as s — > T. This gives rise to the first 
question. 

(Ql)Is there an intrinsic way to relate the growth rate (as s — > T) of the 
derivatives of H to the irregularity of the payoff hi 

The answer will be yes via the notion of fractional smoothness introduced 
below, see Theorems 1 and 2. 



H{t,x)=-E <Q {h{S T )\S t =x) 




(1.1) 
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The estimation of stochastic integrals is usually performed with L 2 -norms, 
but in our financial setting, both measures F and Q can be considered. For 
practitioners, errors under the historical probability P are presumably more 
relevant, while the mathematical treatment under the risk- neutral measure (Q 
is simpler in our context (because the tracking error process (f .1) is a (Q-local 
martingale). 

(Q2) Is the definition of fractional smoothness affected by the choice of a 
specific measure? Do the L 2 - convergence rates depend on the choice of the 
probability measures P or (Q? 

In the context we consider the answer concerning the fractional smoothness 
is usually no in the sense of the comments after Theorem 2. Concerning 
the approximation rates the same is checked for examples so far (see the 
remarks after Theorem 9). 

Beyond the approach to measure tracking errors in L 2l we could alternatively 
identify the weak limit of the re-normalizcd tracking error. 

(Q3) Do the weak convergence rates coincide with those in the L 2 sense? 

The answer is not necessarily, as there are counter-examples in which the 
convergence in L 2 and in distribution hold at different rates, see Section 
1.5. 

Finally, through an efficient choice of re-balancing dates r, one can expect to 
reduce tracking errors and improve the risk management of options. 

(Q4) Which time-nets t = (U)f =0 lead to optimal convergence rates? And 
how to relate them to the fractional smoothness of the payoff? 

As answer we get that according to the index of fractional smoothness 
of the payoff, one can define explicitly re-balancing times achieving the 
optimal convergence rates, see Section 1.5. 

These preliminary questions serve as references for the reader when reading 
the next sections. 

Organization of the paper. 

First, we define the probabilistic framework and the assumptions used through- 
out this work. Then in Section 1.2, we define the fractional smoothness and 
provide basic properties: we choose a presentation that is quite illuminat- 
ing regarding the previous preliminary questions. In Section 1.3, we take an- 
other view on fractional smoothness using the interpolation theory. In Section 
1.4, we consider examples of terminal conditions and identify their fractional 
smoothness. Then, in Section 1.5, we go back to the analysis of discrete time 
hedging errors and state the main results. We close by further developments 
and applications of the fractional smoothness in Section 1.6. 
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Assumptions. 

Let us define the probabilistic setting used in the following. We fix a d- 
dimensional Brownian motion W = (Wt)t£[o,T] defined on a complete prob- 
ability space (/2,.Ft,P) and we let (J 7 t)te[o,T] be the augmentation of the 
natural filtration of W. The log-asset X is the solution of the (i-dimensional 
forward diffusion 

X t = x a + f b(s,X s )ds+ f <j(s,X s )dW s . 
Jo Jo 

To state the results, we mainly consider two types of assumptions: 

(SDE) d > 1 and b, a e C fc °°([0, T] x R d ) and aa* > SI nd for some S > 0. 
(GBM) d = 1 and X t = ln(S t ) =W t - (t/2). 

The smoothness conditions in (SDE) are too strong, they are chosen to sim- 
plify the presentation. Whenever useful to simplify even more, we may con- 
sider the very simple case of the geometric Brownian motion (GBM) (here, 
the asset is a martingale, meaning that P = Q). The reader is referred to the 
corresponding original papers for the possible weaker conditions. 

In the following | • | stands for the Euclidean norm and A ^ c B for A/c < 
B < cA if c > 1 and A, B > 0. Expectations and conditional expectations 
under P are simply denoted by E(.) and E(.|J r t ), while under Q, we indicate 
explicitly the dependency w.r.t. the probability measure by writing Eq(.) and 
E Q (.|^). 

1.2 Definition of fractional smoothness and basic 
properties 

Fractional smoothness on the Wiener space can be defined in various ways, 
see [Wat93, Hir99]. Our approach is motivated by the questions discussed in 
Section 1.1. Since we consider only random variables of the form Z = g(Xx) = 
1i(St) (a function of the process at maturity T), the time T plays a specific 
role in our definition. It would be necessary to modify our definition for more 
general dependencies like Z = g(X tl , ■ ■ ■ ,X tn ), see [GGG10]. 

Definition 1. Assume that Z e L 2 (P). 

(i) For < 6 < 1 we let Z e B 9 2 qo provided that, for all < t < T, 

\\Z-K(Z\F t )\\ L2(P) <c(T-t)-i. 

(ii) For < 9 < 1 we let Z e 2 provided that 

J\t - t)- l - e \\Z - E(Z|^)||| 2(p) rfi < oo. 
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The spaces B^.g above will always be obtained by the conditional expectation 
and the L2-norm under the measure P. Therefore we omit the dependency on 
P in the notation. 

The following properties follow straight from the definition: 

Proposition 1. For < 9 < r/ < 1 and p,q G {2,oo} we have that 

§2,00 B|g and B| 2 CBL. 

Given a bounded 3 measurable g : R d — > R and (t, x) e [0, T) x R d , we let 
u(t,x) :=M(g(X T )\X t = x), 

The following equivalences are useful to exploit properties of Bj 2 and ^2,00- 

Theorem 1 ([GMlOa, Proposition 4]). Under the condition (SDE), for 
< 9 < 1 and a bounded g, the following assertions are equivalent: 

(i) g(X T ) e B^ 2 . 

(ii) r T (T-t)- e E|V x «(t,X t )| 2 di <oo. 
(hi) f o T (T-i) 1 - E|£> 2 M(t,X t )| 2 C ft < 00 • 



Theorem 2 ([GMlOa, Lemma 6]). Under the condition (SDE), for < 
6 < 1 and a bounded g, the following assertions are equivalent: 

{i)g(X T )eB 2iOO . 

(ii) sup^o^T-^ElVxW^Xt)! 2 < 00. 

(hi) For < 9 < 1 we ftaue fAaf sup te[0 T) (T - i) 2 ~ e E |.D 2 w(i, X t ) | 2 < 00. 

Theorems 1 and 2 generalize results obtained in [GG04] and [GH07]. We see 
that the fractional smoothness index 9 measures exactly the growth rate of 
the derivatives of the associated PDE solved by u (see question (Ql) in the 
introduction) . 

The two above theorems are also valid if u is computed using the risk- 
neutral measure Q (i.e. uq(*,x) = Wi^{g{XT)\X t = x)), while the other 
L2-norms are computed under P. For instance, for < 9 < 1 the equivalence 
of (i) and (ii) of Theorem 1 becomes g(Xx) € B^ if and only if J Q T (T — 

3 Here again, the boundedness assumptions on p can be weakened and we refer to 
the original papers. 
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£)~ e Ep \V x u<Q(t, X t )\ 2 dt < oo where we have indicated explicitly if the L 2 - 
norms or conditional expectations are computed under IP or Q. This property 
can be established following [GMlOa] and the proof of [GMlOb, Lemma 7]. 
This accommodates well the fact that the price functions are usually computed 
under the risk-neutral measure, while hedging is made under the historical 
probability (see question (Q2) in the introduction). 

Proof (Simplified proof of Theorem 2). 

We sketch the proof in the simple case where X — W is a linear Brow- 
nian motion, d = 1 and 6 e (0,1). First, (u(t,W t ) = TE,(g(WT)\Tt))t<T is 
a martingale in L,2(P). In addition, for any fixed < S < T the processes 
(V x u(t, W t ))t<T~s and (D 2 u(t, W t )) t <T-6 are i 2 (P)-martingales. This prop- 
erty is obtained by checking that V x u and D 2 u both solve the parabolic heat 
equation and that certain integrability assumptions are satisfied. Then by 
Ito's formula, one obtains for < s < t < T that 

g(W T ) - u{t,W t ) = J V x u(s,W s )dW s , (1.2) 
V x u(t,W t )-V x u(s,W s )= [ D 2 u(r,W r )dW r . (1.3) 

J s 

From the ltd isometry, one deduces from (1.2) that E|(/(M / r) — u(t, Wt)\ 2 — 
J t E|V x w(s, W s )\ 2 ds and it follows that (ii) =>■ (i). Similarly from (1.3) one 
obtains 



E|V*u(t, W t )\ 2 < 2E|V x u(0, W Q )\ 2 + 2 / 1E\D 2 u(r, W r )\ 2 dr 

Jo 

which proves (Hi) => (ii). Finally, we show (i) (Hi). Standard computations 
give that 



., , e 2 < T -*) 



(D " MMr,)= °~ I " : » ^(T-t) dX 



z=Wt 



= / 9(x 



(x - z) 2 - (T -t) e ^=*T 

dx 



r (T - 1) 2 y/2w(T-t) 

(W T - W t ) 2 - (T-t) 
(T-t) 



z=W t 



= E ([ff(Wr) - E(g(^ T )|J- t )] ( ^ T ( y t)2 t)2 (r V * 
which implies that 

||r> 2 «(t,wi)|| L2(P) < H^ 2 | [g(w - r) _ E( g (w T )17-t)|U 3 (p) 

so that we are done. 
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1.3 Connection to real interpolation 

Let us connect Definition 1 to the classical notion of fractional smoothness 
which also explains the notation we have used. Inparticular, this connection 
will make clear the difference between IB^oo an d ^2,2- 

Definition 2 ([BL76, BS88]). Assume a couple of Banach spaces (Eq,E\) 
so that Ei is continuously embedded into Eq. Given ie£o and < A < 00, 
the K-functional is given by 

K(x, A; Eo,Ei) := mf{\\x \\ Ea + AHxiH^ : x = x a + xi}. 

Moreover, given < 9 < 1 and 1 < q < 00 we define the real interpolation 
norm 

\\x\\ e , q := llA-^^A;^,^)!!^^)^) 
and the space (Eo,Ei)g jq := {x G Eq : \\x\\g, q < 00}. 

With our setting (E\ is continuously embedded into E n ) we obtain the fol- 
lowing lexicographical ordering of the interpolation spaces: 

Ei C (Eu,Ei)e, p C (E ,Ei)g^ q C (Eo,Ei) Vir C E n 

for all < 77 < 9 < 1, 1 < p < q < 00 and all 1 < r < 00. 

We apply this concept to the analysis on the Wiener space, which needs to 
introduce some standard notation (see [Nu06, Sections 1.1 and 1.2]). Let H be 
a separable real Hilbcrt space with the scalar product denoted by (., .)h and 
(A4,S, n) be a complete probability space. We assume an isonormal family 
g = {gh ■ h G H} of centered Gaussian random variables, i.e. 

9k) = (h, k) H for all h,ke H, 

and that S is the completed a- field generated by the random variables {gh : 

heH}. 

For each n > 1, we denote by H n the closed linear subspace of 1^2 (m) 
generated by the random variables {H n {gh) : h £ H, ||/i||h = 1} where 

H " w = 1 ^- eT ^ (<rT) ' (1 - 4) 

i.e. the n-th Hcrmite polynomial. "Ho is the set of constants. H n is the so-called 
Wiener chaos of order n and we define by J n : L 2 (/1) — > i 2 (/z) the orthogonal 
projection onto "H„. The following orthogonal decomposition is known as the 
Wiener chaos decomposition: 



n=0 



Now, we are in a position to define the Malliavin Sobolev space and Malliavin 
Besov space. 
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Definition 3. The Malliavin Sobolev space Di^/x) C L 2 {n) is given by 

®M ■■= e L 2 (/x) : II^Id^) := ^f> + 1)|| Jn^||i aM ^ < oo 

Moreover, given < 9 < 1 and 1 < g < oo, we define the Malliavin Besov 
space 

KM ■= (£ 2 (/i),lD li2 (/i))fl,, 
of fractional smoothness 9 with fine parameter q. 

We use this construction in the case that H = t 2 and M. = R d , S is the 
completion of the Borcl er-algebra on R d and /i = 7^ is the d-dimcnsional 
standard Gaussian measure. The family of Gaussian random variables is given 

by 

g h (x) := (x, h) for x e M = JR d and he H = i d 2 . 

To make the connection between the definitions of Bj q and Bj q (jd) for q e 
{2,oo} we let, as before, (W t )te[o,i] be the standard <i-dimensional Brownian 
motion on (J?, J 7 , P, (^t)t£[o,i])- Then we have 

Theorem 3 ([GH07, Corollary 2.3]). For < 9 < 1, 1 < q < 00, and 

g e L 2 (^fd) one has 



£,([0,1),^) 



IMIb» i4 ( 7< i) ~c |MU 2 ( 7d ) + 2 ||Mi -M t || La(P) 

w/iere M t := E (<7(Wi)|.7-i) and c > 1 depends on (9, q) only. 

Applying this theorem to q = 00 gives that 

yWviM ~c IMk( 7d ) + sup (1 - l|Mj - M t || ia(p) , 
whereas q = 2 gives that 

iMk 2(7d) ~c iMU 2(7d) + ii m i - M *nL(P) dt ) 

which brings us back to Definition 1. 



i=l, 
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Multi- dimensional Black- Scholes- Samuels on model. 

This is a log-normal model which dynamics on the price and the log-price can 
be written as 

dS\ = SI ^ o-ijdWi + mdt^j , l<i<d, 
d 

X\ = log(4) +^<r#/ + ( M i - - O o\)t, 

where a t :— ^Jj2j a ij ■ Assume that (<7ij)f - =l is invertible. To the payoff func- 
tion S i y h(S) 7 we associate g(x\, Xd) :— h ^ ^s^e^ 1 ^ 1 a ^ X3+fl ' 
From this we see that 

J£<(7J) if and only if h(5i) G B^, 
for q e {2, oo} and g G L 2 ( r )d)- 
Remark 1. In the case 0=1 we get that 

<? € E>i, 2 (7d) if and only if e b£ i00 

for all g G L 2 (jd)- This can be checked by using arguments from the proof of 
[GH07^ Corollary 2.3]. 



1.4 Examples 

In this section, we provide examples of random variables Z — g(Xx) for which 
we determine the fractional smoothness. 

Example 1 (Lipschitz function). The case, where the fractional smoothness is 
obvious, is the Lipschitz case. Assume a Lipschitz function g : R d — > R with 
constant L > 0, i.e. \g(x) — g(y)\ < L\x — y\ and assume (SDE). Then one 
has that 

E \g(X T ) - TS(g(X T )\T t )\ 2 < E \g(X T ) - g(X t )\ 2 

< L 2 ~E\X T -X t \ 2 

< L 2 c 2 (T-t), 

using standard estimates on the increments of X. Hence, g(Xx) G B2 00 ■ This 
example includes call and put payoffs, i.e. g(x) = (x—K) + or g{x) = (K—x) + . 
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Exactly the same argument as above yields for 0-Holdcr functions g with 
9 e (0, 1) that g{X T ) e B^.oo- But the situation is here not as clear as one 
expects as shown by 



Example 2. Assume the setting (GBM) and that 

h e (x) := (x-K) e + 



for some K > and < < 1/2. Then it is shown in [GT01, Lemma 2] 
(under more general assumptions) that E|D 2 w(i, X t )\ 2 < c(T — t)- 3 / 2 + e S o 
that Theorem 2 gives that 

For 1/2 < 9 < 1 one gets h e {S T ) € B£ i00 . 

Example 3 (Binary option). Generally, indicator functions yield to a fractional 
smoothness of order i. In the case X — W, d = 1 and g(x) = 1[l,oo)( x ) with 
LgE one has 



and <7(Wt) <E B| ^ because of Theorem 2. This can be extended to the (SDE) 
follows: Our assumption guarantees that X has a transition density r 
such that 



for some k > and all < s < t < T, where j t is the Gaussian density with 
zero expectation and variance t (see [Fr64]). Then we can compute that 



E |1[l,oo)(*t) - E (1 [l ,oo)(^t)|^) I 

< E |l[ Li00 )(X T ) - l[ Li00 )(X t )| 

= P(X T < L < X t ) + P(X t < L < Xt) 

< k, 2 [F(W kT <L-x < W Kt ) + V{W Kt <L-x < W kT )\ 

< cxjT-t 



where Xq = xq so that 1[£ i(X) \(.Xt) € IBloo- The application in financial 
mathematics is done via S t — e Xt which gives, for a positive strike K > 0, 




so that 



nv xU (t,w t )\ 2 



i 
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1{S T >K} - l{X T >1o g K} € B| iOCJ . 

In our context the fractional smoothness of jump functions (under differ- 
ent assumptions) was considered in [GT01, Gei02, GG04]. In certain multi- 
dimensional settings one can deduce for g(x) = i{ Xl >K 1 .---,x d >K d } ( or variants 

of it) the same fractional smoothness from the 1-dimcnsional case. Finally, the 

~ i 

indicator function g{x) = 1_d(x) of a C 2 -domain D also leads to g(Xr) € oo 
(see [GM05, Proposition 1.2]). 

Example 4 (an extreme case). By the choice of the previous examples, we 
emphasize that random variables g{Xr) — /i(St), usually used in financial 
applications, belong to a space B2 x for some 9 € (0, 1]. However, it is not 

true that U ee ( 0i i]lB2 iOO = ^(P). The following result gives a way to con- 
struct g(Wi) belonging to ^2(P) (here W is the linear Brownian motion) but 
5(^1)^^2,00 for all 0e (0,1]: 

Proposition 2 ([GH07]). Let < 6 < 1, g = J2T=o a k H k e £2(71), where 
(Hk)k>o is the orthogonal basis of Hermite polynomials defined in (1.4). Then 

00 

g(Wi) e B200 */ and only if sup (1 - t) 1 " 8 V kt k ~ x a\ < 00. 

o<t<i k=1 

Approximation properties as described in Section 1.5.2 for g with g(Wi) <G 
^2(P) \ Uo<e<i ^2.00 w ere studied in [Huj06] and [Sep08]. 



1.5 Applications 

In this section we discuss some applications in stochastic finance which lead 
us to the fractional smoothness as introduced above. As mentioned at the 
beginning, a central role is played by the tracking error that arises when 
discrete time hedging is used, instead of a continuous time strategy. For the 
sake of convenience, we briefly recall the notation: 

• the option payoff at maturity T is Z = h(Sr); 

• the fair price function is H(t,x) = M<Q(h(ST)\St = %)', 

• the n re-balancing dates are defined by a deterministic time-net r = (U)^ =0 
with = t Q < h < ■ ■ ■ < t n -i <t n = T; 

• the resulting tracking error process C{Z 1 t) = (C t (Z, r)) tG [ |T j is given by 

n-l 

C t (Z, t) := E Q (Z| T t ) - E Q Z - ]T V x H{t u S u ) ■ (5 ti+lAt - S UAt ). 

i=0 
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1.5.1 Weak limits of error processes 

Weak limits of stochastic processes have been intensively studied in the 
literature; see, for instance, [KP91, Jac97, JS03]. For the particular prob- 
lem of the weak convergence of the tracking error the reader is_ referred to 
[R008O, GT01, HM05, GT09]. To formulate our results, we let W = (W t )t>o 
be a standard Brownian motion starting at zero defined on some auxiliary 
probability space, where we may and do assume that all paths are continu- 
ous. In the following =^>c[o.s] stands for the weak convergence in C[0,s] for 
some s > 0. 

In this paragraph we assume that T = 1 4 and that S is the standard 
geometric Brownian motion, i.e. the setting of (GBM) and P = Q. The 
following result is the starting point of this section: 

Theorem 4 ([GT01]). Let r n = (i/n)™ =0 be the equidistant time-nets and 
let Z := 1[k,oc)(Si) be the payoff of a digital option with strike price K > 0. 
Then one has that 



dt 



y/nCi(Z,T n ) =^> W 1 rU q2a 2 H 

3 Jo | 5 t^(*' S * 

where => denotes the weak convergence as n goes to infinity. 



The remarkable fact is that the weak limit is not square-integrable. In the 
following we describe a way to increase the integrability of the weak limit. This 
is of particular interest for risk management purposes, as a higher integrability 
gives better tail-estimates. The idea is to use adapted time-nets that are more 
concentrated close to maturity. They are defined as follows: Given a parameter 
9 e (0, 1], we define the nets r" ,e by 

k 
n 



C := 1 - I 1 - - 



For 6 = 1 we have the equidistant time- nets, i.e. i£' = K Now we have 



Theorem 5 ([GT09]). Let < 6 < 1, Z = h{S 1 ) € L 2 (P) and < s < 1. 
Then 



t€[0,s] 



(V^C t (Z,T n ' )) t6[M =>C[0,,] \s*f£(r,S r )\'dr 

Moreover, the following assertions are equivalent: 



4 With T — 1 we are in accordance with the quoted literature that used Hermite 
polynomials. Of course, we could do a re-scaling to T > afterwards. 
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(i) One has ft(Si) G B e 2 2 for < 6 < 1 or ft(Si) G B^ )00 for 6 = 1. 

(ii) On some stochastic basis there exists a continuous square-integrable mar- 
tingale M = (M t ) te [ 0l i] such that y/nC(Z,T n - e ) =>c[o,i] M - 

(iii) For 

A -J Q 29 bt dx* {t ' St) M 
one has that E^4 < oo and 

The theorem above gives us one way to consider the L p -setting for 2 < 
p < oo. Given a diffcrentiable function ?/> : (0, oo) — >• R we let 

(A^)(a;) := xtp'(x) - ip(x). 

In the following AH(t,x) means that A acts on the a;- variable of the function 
H(t,x). 

Definition 4. For h(Si) G L 2 (P), < < 1, and < t < 1 we let 

Df^hiSt) := ^ /"'(I - u)- 1 ^ [AH(u A i, S uAt ) - AH(0, S )]du. 
1 Jo 

For 6 = 1 and i G [0, 1) we let D? !l ft(Si) := Aff (i, S t ) - Aff(0, S ). 

The process D S ' 6 h(Si) = (Dfh(Si)) te [ 0} i) is a quadratic integrable mar- 
tingale on the half open time interval [0, 1). Using the Riemann-Liouville op- 
erator of partial integration the process D S ' 6 h(Si) can be interpreted as a 
fractional differentiation of order 8 in x (sec [GT09]). The point of the con- 
struction of D S ' 9 h(S\) is that we may have L p -singularities of St §§■(£, S t ) as 
t 1 1 whereas D S ' e h(S 1 ) remains -Lp-bounded. 

Theorem 6 ([GT09]). For 2 < p < oo, < 6 < I, and Z = ft (Si) G L 2 (P) 
the following assertions are equivalent: 

(i) On some stochastic basis there exists a continuous L p (JP) -integrable mar- 
tingale M such that y/nC(Z,T n ' 6 ) =^c[o,i] M - 

(ii) The martingale D S ' 8 h(Si) is bounded in L p (P). 
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1.5.2 Z/ 2 -estimates of the tracking error 

In this section we work in the 1-dimcnsional martingale case assuming (SDE) 
with cr(t, x) = a{x) and b(t, x) — — \a 2 {x) (meaning P = Q). The payoff 
function h is polynomially bounded and the option maturity is T > 0. We 
remind the reader about the time-nets T n ' S given by 

C-r(i-(i-|)*) 

and that for 9 — 1 we obtain the equidistant nets. Let us first check what 
quadratic hedging error one can expect at all if the portfolio is re-balanced 
n-times. The answer is the rate 1/y/n as shown by 

Theorem 7 ([GG04, Theorem 2.5]). Assume that there are no constants 
co, ci e It such that ft (St) = co + c\St a.s. Then 

Jnf ni\\C(h(S T ), (t fc )JLo)lk(P) > 

a=t <-<t n = T 

where the inftmum is taken over deterministic time-nets. 

This was extended to the case of random time-nets in [GG06] in the case 
of the geometric Brownian motion. 

Now we continue with the case of equidistant time-nets which are often 
used in discretizations. 

Equidistant time-nets. 

Here a starting point is the following result of Zhang: 

Theorem 8 ([Zha99, Theorem 2.4.1]). Assume that h : R -> R is a 
Lipschitz function. Then we have that 

limn^||C(ft(S T ),T^ 1 )|| i2(P) e[0,^). 

This is the result one would expect: Given a Lipschitz payoff, the Li-x&te, of 
the error is 1/2 for equidistant nets. But this is not the case in general as 
shown in 

Theorem 9 ([GT01, Theorem 1]). For h(x) = l [Ki00 )(x) for some K > 
we have that 

limn^C^T),^' 1 )!!^) e (0,oo). 
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This means that the L 2 -approximation rate for the binary option is n 1 / 4 if 
one uses equidistant nets. The two above results also hold true for appropriate 
Qt^P (i.e. S is not martingale) where the outer L 2 -norm is computed w.r.t. 
the historical probability P (cf. the remarks after Theorem 2). 

Theorems 8 and 9 lead naturally to two questions: What is the reason for 
the rate 1/4 and, secondly, can one improve the rate 1/4? Both questions can 
be answered by the usage of the concept of fractional smoothness. 



Theorem 10 ([GG04, Theorems 2.3 and 2.8]). For < 9 < 1 and a 

polynomially bounded h : (0, oo) — >• 1R the following assertions are equivalent: 

(i) h(S T ) e B% t00 . 

(ii) sup n n«||C(/i(Sr),r n ' 1 )|U a( p) < oo. 

In particular, it turns out that h(Sr) <E E>i ;2 if and only if 
supniHCWSTj.T"- 1 )!!^) <oo, 

n 

see [GG04, Theorem 2.6], where ©1,2 is the Malliavin Sobolev space obtained 
from the construction in Section 1.3 with H = L 2 [0, T] and gh := Jq h(t)dW t . 

For the binary option one has in Theorem 10 that 9 = 1/2 (cf. Example 3 
in Section 1.4). This recovers the rate 1/4 obtained in Theorem 9. 

Non equidistant time-nets. 

Next we show how to obtain the optimal rate n 1 / 2 by a suitable choice of 
the trading dates (see question (Q4) in Section 1.1). We can combine [GG04, 
Lemmas 3.2 and 5.3] and [GH07, Lemma 3.8] to get 

Theorem 11. For < 9 < 1 and a polynomially bounded h : (0, 00) — > R the 
following assertions are equivalent: 

(i) tf(T-ty- e v\s?gt(t,s t )\ 2 dt<^. 

(ii) sup„n^||C(MS T ),T"< e )|| i2(P) < 00. 

For < 9 < 1 (and at least a bounded h) the condition Theorem ll(i) is 
equivalent to 

(i') h(S T ) e B e 2 . 2 

which can be checked by using Theorem 1. For the binary option this gives 
that 

supni||C(l [Ki00) (S*T),T™< r ')|| L2 ( P ) < 00 
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for any strike K > and < rj < 1/2. 

For the next two theorems we assume that T = 1, that S t = e w *~i and 
that h might be general, i.e. not polynomially bounded. The formulation of 
Theorem 1 1 in the language of the interpolation spaces introduced in Section 
1.3 gives 

Theorem 12 ([GH07, Theorem 3.2]). For < 6 < 1 and G L 2 (F) 

the following assertions are equivalent: 

(i) fc(e'~ (1/2) ) € Bg i2 (7i) i/0<6><1 anrf/i(e-( 1 /2)) eE)li2 ( 7l ) #0 = 1. 

(ii) sup n r"' e )|| i2(P) <oo. 

And Theorem 10 can be extended in this context to the full scale of real 
interpolation spaces as 

Theorem 13 ([GH07, Theorem 3.5]). For 1 < q < oo, < 6 < 1 and 

h(Si) € L 2 {P) the following assertions are equivalent: 

(i) h(e-VM) € B^(7i). 

/ 6 1 \ ^ 

(ii) n5-iaj < oo for a n := ||C(/i(S'i),T ,l4 )|| i2(P) . 

Concluding remarks 

(i) The higher dimensional case for X was considered in the literature as 
well. Roughly speaking, one can analogously obtain upper bounds, how- 
ever precise lower bounds as in the one-dimensional case are still missing. 
This is due to the fact that a characterization of the L 2 -error proved in 
[Gci02, Theorem 4.4] and [GG04, Lemma 3.2] is missing for higher di- 
mensions. However, after Zhang [Zha99] started with the regular case, 
Temam [Tem03] extended results from [GT01] to higher dimensions and 
Hujo [Huj05] used non-uniform time-nets to improve the approximation 
rates for certain irregular payoffs to the optimal rate l/^/n in this set- 
ting. 

(ii) Seppala [Sep08] found a criterion to characterize under certain conditions 
that there is a constant c > such that 

r=( mf o ||C(MSi),r)|U 2(P) <-^ 

o=t <-<T n =i 

where deterministic time-nets are taken. It should be noted that one has 
a non-linear approximation problem as the time-nets may change for 
fixed n from payoff to payoff h. 

(iii) In the above discussion, the time-nets r are deterministic. Alternatively, 
one can allow the time-nets to be stochastic and adapted. This issue has 
been handled by [MP99] using optimal stopping tools. The estimation of 
convergence rates is an open question. However, it was shown in [GG06] 
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that the random time-nets do not improve the best possible approxima- 
tion rate 1/y/n in the case (GBM) when in the n-th approximation a 
sequence of n stopping times is used, 
(iv) Similar studies can be performed when studying the Delta-Gamma hedg- 
ing strategies. Instead of hedging the payoff using only the asset, we use 
other traded options written on the same asset. For a one-dimensional 
asset, if the price of the additional option is (P(t, S t ))o<t<T, the num- 
bers of options P and assets to hold at time U are respectively equal 
to 

P d 2 s H{ti,S ti ) s djH(ti,S ti ) . 

5t * - d>p(u,s u ) and 5 ** - dsHt *> St *) ~ a*p(u,s u ) dsPiu ' Stt) - 

In [GMlOb, Theorem 6], considering a multi-dimensional Black-Scholcs 
model, it is established that for an exponentially bounded payoff such 
that g(Xr) € ^2oo for some < 6 < 1, the use of equidistant time- 
nets leads to the same convergence rate l/n 9 / 2 as for the delta hedging 
strategy. On the contrary, the use of non equidistant time-nets T n ' v with 
< T) < 9/2 enables us to obtain the improved convergence rate l/n. 



1.6 Further developments 

1.6.1 Backward stochastic differential equations 

Makhlouf and the second author applied in [GMlOa] the concept of fractional 
smoothness to backward stochastic differential equations of the type 

Y t = g(X T ) + £ f(s, X s ,Y s ,Z s )ds - J* Z s dW s 

where X = (X t )te[o,T] is our forward diffusion and the generator / is continu- 
ous in its four arguments, continuously differentiable in (x, y, z) with uniformly 
bounded derivatives. These equations are particularly useful in stochastic fi- 
nance, since they allow to take into account market frictions and constraints 
(we refer to [EPQ97] for a more complete account on this subject). 

Solving numerically this type of equation is a challenging issue since it 
concerns a non- linear problem (due to the generator /), generally defined in 
a multi-dimensional setting. One possible approach consists in approximating 
the BSDE using a discrete-time dynamic programming equation (see [Zha04, 
BT04, LGW06] among others). One of the main error contribution is related 
to the L 2 -regularity on Z, defined by 

£(Z,t)=J2 WZt-Zt^Wl^dt. 
»=i Jt «-i 
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If / were equal to 0, then the Z-component is given by z t = V x u(t, X t )a(t, X t ) 
where u(t,x) = TE(g(XT)\X t = x). Studying the Z^-regularity of z is thus 
very similar to the analysis of the tracking error presented in Section 1.5. 
Additionally, using BSDE techniques, one can prove explicit upper bounds 
for the difference Z — z 



Theorem 14 ([GMlOa, Corollary 14]). Assume (SDE) and g(X T ) E 
Bg.oo for < 6 < 1. Then, for some c > 0, one has that 



z t \ < c 



E \g(X T )-lE(g(X T )\T s )\ 2 \F t 



T-s 



■ds + c{T-t), 



M\Z t -z t \ 2 < c(T-t) 6 



Taking advantage of this approximation result close to the time singularity, 
we can prove that the estimate of £(z,t) (linear case) transfers to £{Z 1 t) 
(non-linear case) and get 



Theorem 15 ([GMlOa, Theorem 21]). Assume (SDE), g{X T ) e 
that < i] < 9 < 1 or n = 9 = 1 . Then one has that 



£{Z,T n ") < -. 

n 



In [GGG10], extensions of the above in different directions are discussed. 

1.6.2 Levy processes 

An extension of the results of [GH07] to Levy Processes is done by C. Geiss 
and Laukkarinen in [GL10]. Moreover, Tankov and Broden proved in [TB09] 
results along the line of [GT01]. 

1.6.3 Multigrid Monte-Carlo Methods 

In the context of Multigrid Monte-Carlo Methods it turned out that the con- 
cept of fractional smoothness is useful as well. The reader is referred to the 
papers of Avikainen [Av09a, Av09b]. 
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